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Abstract 



We address the problem of solution of the QCD three-particle evolution equations which govern 
the Q 2 -dependence of the chiral-even quark-gluon-quark and three-gluon correlators contributing 

"xl" ' 

t~"-> I to a number of asymmetries at leading order and the transversely polarized structure function 

O '• 

Q\ ■ ^(^Bj)- The quark-gluon-quark case is completely integrable in multicolour limit and corresponds 



to a spin chain with non-periodic boundary conditions, while the pure gluonic sector contains, 



Ph| apart from a piece in the Hamiltonian equivalent to XXX Heisenberg magnet of spin s = — ~, a 



1 

Ph. non-integrable addendum which can be treated perturbatively for a bulk of the spectrum except 
for a few lowest energy levels. We construct a quasiclassical expansion with respect to the total 
. £h ! conformal spin of the system and describe fairly well the energy spectra of quark-gluon-quark and 
^ . three-gluon systems. 
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1 Twist-three effects in high energy scattering. 



The naive parton model [1] provided a first successful microscopic description of several high energy 
inclusive processes and endowed them with intuitive probabilistic interpretation as the scattering 
of a probe on an incoherent bunch of free and collinear constituents (quark and gluons) in a hadron 
target by expressing the cross sections in terms of the densities of partons in a (incoming) hadron 
or hadrons in a (outgoing) parton. Later it was realized that this is only the first asymptotic term 
in the infinite series in inverse powers of a hard momentum scale Q where all quantum mechanical 
interference effects enter the game as power corrections and, loosely speaking, can be interpreted 
as an overlap of hadron wave functions with different number of partons. For the deep inelastic 
scattering this series has a firm field theoretical ground in the form of the operator product 
expansion (OPE) || [|. For other reactions which do not admit the OPE a set of factorization 
theorems |3J] has been proven within the framework of QCD, however, not far beyond the first non- 
leading term |5| which goes under the name of higher twist. These are the latter terms which are 
of particular interest since being correlations of more than two field operators they can give new 
insights into the QCD dynamics in the non-perturbative domain which is still out of a systematic 
theoretical control. 

Typical higher twists in hard reactions are usually associated with power suppressed contribu- 
tions in a large momentum scale on the background of a dominating twist-two cross section and 
while being important cannot be studied experimentally with a high accuracy. However, it is not 
always the case as there exists a class of transverse spin dependent processes where the relevant 
observable is an asymmetry, Act = <r(s±) — cr(— s±), (not a cross section a = cr(s±) + c(—s±)) 
and where the twist-three contributes as a leading effect not contaminated by a lower twist. The 
well-known twist-three structure function ^(^Bj) H measured in the polarized deep inelastic scat- 
tering is their good representative. It acquires an operator definition via a Fourier transform of 
the matrix element of a non-local operator of the quark fields separated by a light-like distance?] 

1 f d k — 

g T (x) = 9l (x)+g 2 (x) = - / e^l^Oh^MI/*)- (1) 

2 J 271 

Being of the form similar to the conventional unpolarized distributions this quantity is implicitly 
interaction dependent. This can be made manifest by virtue of QCD equations of motion and the 
Lorentz invariance. Its genuine twist-three part can be reduced in leading order in the QCD cou- 
pling constant to an integral of a quark-gluon-quark correlation function (whose precise definition 

2 We have omitted here (and everywhere below) the path ordered exponential assuming the light-cone gauge 
B+ = 0. 
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is given in Eq. ( |139| )) depending on two momentum fractions Y(x,x') ~ (i/jGi/j), namely 

doc P doc 



(2) 



Evaluating the correction in a s to the structure function g2 requires the introduction of three- 
gluon operators GGG (see Eq. (||)) which enter through the quark loop couplingP] [[| and might 



be responsible for the small-a^j behaviour of ^(^Bj), — the issue which is under debate |TT 



Another set of observables is single transverse spin asymmetries in hadronic reactions where 
only one particle in the initial state is polarized (or only polarization of a single particle in the final 
state is tagged), i.e. the inclusive Compton scattering •yp — > 7X |12|, the direct photon |13|, fUjl , 
pion |13| or jet production in pj>-collisions, the Drell-Yan process with measured azimuthal angular 
dependence of a lepton |15| etc. They are of great interest since they involve almost the same 
quark-gluon-quark and three-gluon correlation functions discussed above but not folded with a 
coefficient function and thus depend on both (although equal) momentum fractions. This can 
be used for their extraction or at least for a direct confrontation with phenomenological models 
available. 

Since the data are normally can not be taken at the same values of the hard momentum variable 
Q, the correlation functions which enter cross sections are measured with different resolutions and, 
thus, they differ for different scales. Therefore, the question naturally appears about their relation. 
Their dependence on Q is only logarithmic and within the context of QCD it results from the 
behaviour of the theory at the cut-off. This can be consistently described within the framework of 
the renormalization group for composite operators. In the language of the QCD improved parton 
model the task of its description is translated into the construction and solution of the evolution 
equation for the multi-parton correlation functions. At leading order in the coupling constant, 
to which we restrict ourselves, the former problem is relatively simple and can be resolved in a 
straightforward manner [TjJ [17], [Tj|. However, the second issue is by far more complicated since 
we face here a three-particle problem and a priori there is a little hope to find analytical solutions 
to it. The ultimate result of the scale dependence of a correlator is expected to be of the form 



F(x J x'\Q 2 ) = J2*{ a }(^x') ( 



as(Ql) \ fic)£{a}/Po 



((F {a} (Qt))h (3) 



where £{ a y are the eigenvalues of evolution kernels (and /( c ) is an extracted colour factor) and 

<^{ Q }(x, x') are the corresponding eigenfunctions parametrized by a set quantum numbers {a}. As 

3 Notc, that the results of this reference are not complete since the contribution from twist-three two-gluon 
operators which are related to the three-gluon ones |ic|l were not taken into account and must be included for 
complete treatment. 
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usual Po = ^TpNf — yC*A is the leading term of the QCD /3-function and {{^F{a}{Ql))) stand for 
reduced matrix elements of local operators at a low normalization point Qq. 

Very recently a guiding principle has been found which has allowed to cope the second problem 
and it is based on the integrability of spin lattice models to which the corresponding evolution 
equations can be reduced (19], |22|. These findings are analogous to a similar equivalence 

encountered in the problem of description of the Regge behaviour of QCD amplitudes [^3|, |24|, ^5], 
[26|j . The studies have been undertaken for the three-quark [|l^, |2(| and chiral-odd quark-gluon- 
quark [19|, ^1], 22] operators and allowed to find the anomalous dimensions, hereafter referred to as 
energies, of the eigenstates as well as their eigenfunctions in a WKB type manner where the role 
of the Planck constant has been played by the inverse total conformal spin of the system. These 
results compared quite well with an explicit numerical diagonalization of the mixing matrix and, 
therefore, they provided a reasonable approximation. 

In the present paper we continue the study along this line and give almost complete description 
of the three-gluon and quark-gluon-quark chiral even sector. Our presentation will be organized 
as follows. Next section is devoted to the study of the purely gluonic evolution equation. This 
means that we discard completely from our consideration the mixing with operators containing 
quark fields. We construct the evolution equation for corresponding correlator and because of its 
complexity we are forced to find an effective approximation to it. By making use of the conformal 
symmetry the problem is reduced finally to a quantum mechanical problem for particles with a 
conformal invariant pair-wise interaction and it turns out that the corresponding total Hamiltonian 
is a sum of an integrable piece equivalent to the Heisenberg spin chain of spin s — — | and 
an addendum which breaks the integrability but still can be treated as a small perturbation 
for a large part of the spectrum. However, it plays its crucial role for the generation of a fine 
structure of lowest energy levels. Next, we address the non-singlet quark-gluon-quark sector which 
in the limit of infinite number of colours is shown to be identical to an inhomogeneous spin chain 
with non-periodic boundary conditions. In both cases, gluon and quark-gluon problems, we are 
entirely interested in the polynomial solutions of the lattice models, which play an exceptional and 
distinguished role for physics being implied by the OPE, since each polynomial will correspond to 
a multiplicatively renormalizable local operator. However, non-polynomial solutions are equally 
interesting albeit their physical relevance within the present context remains obscure presently. 
Finally, we summarize. 
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2 Gluonic sector. 



In this section we address the question of the diagonalization of the purely gluonic twist-three 
evolution equation for a correlation function which contributes to g% . We will do it in parallel 
with an operator T which possesses exactly solvable interaction and thus is very insightful for the 
realistic case. Thus, we introduce^] 

^ ) Ka, «s) = 8 I fahc9 ^ \ G^{K,n)G h ^{K 2 n)Gt p {K in ), (4) 

where the totally symmetric, w.r.t. the independent permutation of {a, (3, 7} and {//, z/, p} indices, 

tensor isQ r a/ 3 7;M1/p = V(a, (3, i)g^ a {g^g^ + g^gjp) - \V{ji, v, p)g^ v {g^g p \ i + g^gjp + g^gja)- The 

momentum fraction space functions are 

F(x 1 ,x 3 ) = J ^^ e ^-^{h\H^^)\h), (5) 

with T = Q,T and possess the symmetry properties G(xi,x 3 ) = —G(—x 3 ,—Xi), T(xi,x 3 ) = 
T(— s 3 , —x\). They obey the evolution equation 

d F( Xl ,x 3 ) = -^J 11*45(4 -4 + 4)^ F ({^}IKJ)^ , i,4), (6) 



rflnQ 



where X2 = X3 — xi. 

Since the existing studies || did not provide an insight into the structure of the corre- 
sponding kernel K F we reanalyze the issue anew showing a simple structure of the result which 
allows for a reduction of the problem to a lattice model characterized by a Hamiltonian which 
consists of two parts: exactly integrable piece equivalent to the generalized Heisenberg spin chain 
and a one which violates this property and is responsible for the formation of a fine structure of 
low lying levels. 

2.1 Three-gluon evolution equation. 

To start with, note first of all that only 'good' transverse components of the fields enter the 
operators and, therefore, the latter belong to a special class of the so-called quasi-partonic 



operators [16]. They are distinguished by the property that the renormalization does not move 



4 Thc /abc-coupling in T-operator gives identical zero provided the total derivatives are irrelevant like for the 
forward matrix elements considered presently. Obviously, due to negative C-parity the d a (, c -coupling does not 
contribute to the deep inelastic scattering cross section with an electromagnetic probe. 

5 Here the action of the operator of cyclic permutation V is defined as V(l, 2, 3)/(l, 2, 3) = /(l, 2, 3) -1-/(2, 3, 1) + 

/(3,1,2). 
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Figure 1: Diagrams (in the light-cone gauge) giving rise to the gluon-gluon two-particle evolution 
kernel. The graph with the crossed gluon line corresponds to the contact-type contribution arising 
from the use of the equation of motion. 

them outside of the class and that at leading order the total evolution kernel is a linear combination 
of the twist-two non-forward kernels in subchannels. 

The pair-wise gluon evolution kernel can be decomposed as follows into Lorentz tensors 

"Ka'V;n'v' { X U X 2|^l, X 2 ) 

= ~ {^9iv9i>v>"K V + \e L ^,J°K A + t^t^.^K^ (xx, x 2 \x[, x' 2 ) , (7) 

with g^ u = g^ - n^nl - n*n u , e^ v = e Mi , pCT n*n CT and T^. pa = \ (g^ p gt + gj^g^ - 9^9 p*)', and 
colour structures 

{ 99Kl Y a b , b , = jjrZi6*8*v ag K i o.) + Y c f*°f*v° ggKi M + J^d ah J a , b ,J>3K\, s) + -.., (8) 

where the dots stand for the decuplet and 27-plet contributions which are irrelevant for our 
consequent discussion. 

We can easily evaluate from diagrams in Fig. |l] the kernels for the parity even and transversity 
channels 



1^2 



f,}"Kj c) (x 1 ,x 2 \x' 1 ,x' 2 ) = k(x 1 ,x 2 \x' 1 ,x 2 ) + 2 XlX2 / p^- Q[ n (x 1 ,-x 2 ,x 1 - x[) 

f XiX 2 X\X' X + X 2 X 

x[x 2 (xi + x 2 ) 2 



X\X 2 X\X 1 -\- X 2 X 2 q \ i / / /\ / n \ 



/( C ) " K {c) {xi,x 2 \x' 1 ,x' 2 ) = k(x l ,x 2 \x' 1 ,x' 2 )±(x' 1 ^x 2 ), (10) 

with "+" sign for /( 8s ) = /(x)/2 = and "— " sign for /( 8a ) = Ca/4. For the parity odd sector 

we have 

f{c)" K t) (^I'^ki, x' 2 ) = k(x 1 ,x 2 \x' 1 ,x 2 ) + 2 XlX2 ^^ e^fxi, -x 2 ,xt - x\) 

X\X 2 

+ 2 X flQ u (x 1 ,-x 2 )T(x > 1 ~x' 2 ), 

X\X 2 
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m.c. 



Figure 2: Structure of the total evolution kernel for the three-gluon correlation functions 
G(xi,x 3 ). The symbol m.c. stands for the mirror symmetrical contributions of the last two 
diagrams. 



with "±" signs assigned vice versa to the previous case. The generalized step functions look as 
follows 



Q u (x 1 ,x 2 ) 
Qi 11 (x 1 ,x 2 ,x 3 ) 



e{xt) 9{x 2 ) 



X\ — x 2 x 2 — X\ 

X\d{x\) 



+ 



x 2 9{x 2 ) 



+ 



x 3 9(x 3 ) 



(x 2 - x 1 )(x 1 - x 3 ) (xi - x 2 ){x 2 -x 3 ) (x 1 -x 3 )(x 3 -x 2 )' 
We have introduced the notation 



Xl 
X i 



X\ — X\ 



-yQl^X^Xl - X[) 



+ 



x 2 

X 2 



X 2 



X 2 — X 2 



T Q° u (x 2 , x 2 - x 2 ) 



1 (3o 



for the part which appears in all kernels. 
Introducing the conventions 



:i2i 



and iT i±: ' = | (fO ± iP), we can contract the generalized kernels (|7|) with the tensor structure of 
the operators in question (^) and construct easily the total evolution equation for the three-gluon 
quasi-partonic correlator G (see Fig. [2] for its pictorial representation) 



K 



G 



2™K£+f(x 1 ,x 2 \x[, x' 2 ) +2^Kf+J(x 2 , -x 3 \x' 2 , -x' 3 ) 
+ 99 KV A) (x u -x 3 \x[, -x' 3 ) - ^Kl^{x X: x 2 \x^ -4) - «K { 
- \p S(x 1 - x'^S^ - x' 3 ), 



V—T( I / I \ 

A ) {x 2 , —x 3 \x ± , —x 3 ) 



(13) 



and for the T 



K T = ^ s) (x 1 ,x 2 \x[,x' 2 )+^ s) (x 2 ,-x 3 \x' 2 ,-x' 3 ) + 
- \(3 5{xi - x\)5{x 3 - x 3 ). 



(14) 
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Defining the moments as 

* - 1 iM&n*,*.), (is) 



with the properties Gj = — (—l) J Gj_j, Tj 7 = (— l) J Tj_j we can write, given the results of 
Eqs. flI3p and ([TJ|), the evolution equation for the moments, i.e. for the local operators Fj 7 oc 
gG + ±(id + ) j G + ±(id + ) J ~ j G + ±, in the form: 

F? = -^zY,rUJ)F k J , rf k (J) = C A yf k (J) + (3 8 jk , (16) 

k=0 



d\nQ 2 1 4vr ^ 3 



with the anomalous dimension matrices 

1%{J) = 6 jk {mj + 3) - 3^(1) + ((J + l)C 3 J+ 2 + 3(j + 2)C} +1 ) - 

(j + l)(J-j + l) _ 1 + (-!)' ] 
(J + 2) 3 (j + 2) 2 J 

w r l j eg; {(j+w+2) + (j+2)(j- J + i) (j+2) 3 

(fc + l)(J-J + l) 1 ^ 2 J-fc + l / 



(J + 2) 3 J + 2^ (J- J + l) 3 V ^ 7 '-i+ 2 

+2-t^ ((j + 1) (C* +2 - Cj + (j + 2) (3C* +1 - V + ±)c-]z) 



-2 



(J i~^ 1)4 {{J ~ J + + 3( J ~ J + 2)C*H) 

1 ^i)^-^^l}-fr-!--il 



anc 



7lW- ^{3^(i + 3) -3^(1)} 



- dw — 3 \ { - l) s + (j- J+ i) 2 r^j-j' (18) 

where we introduced the Pochhammer symbol (j)^ = F p^"^ ; the binomial coefficients = 

r(j-H)r{j-j+i) anc ^ ^ e s ^ e P f unc ti° ns — {1> if J > k; 0, if j < A;}. Our mixing matrix 
( |TTD differs from the one obtained in Ref. |8| since the evolution equation was written there for 
a different quantity, i.e. a linear combination of G-functions considered here (|4|,§). But both 
matrices have indeed identical eigenvalues with different eigenvectors. 

6 To get these anomalous dimensions we assumed the gluons being different so that there is no symmetrization 
in the pair- wise kernels. Identity of gluons will just pick up eigenstates with a definite symmetry. With the matrix 



(18) at hand we have the states with all possible symmetries. 
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2.2 Reduction to ID lattice model. 

The key observation is that not all of the parts of the anomalous dimension matrix fllTD play 
an equally important role in the generation of the energy spectrum. Note first that due to the 
antisymmetry of the three-particle correlator G(x\,x 3 ), w.r.t. the interchange of the momentum 
fractions of the first and last gluon, the leading rightmost 4-singularity in the j-plane of the 
vector kernel disappears. For the formation of the upper part of the spectrum j,k ~ J only 
— x'j)] + -distributions are relevant while the fine structure of low levels is generated by the 
K A + i^ T -part of the total kernel. Making use of these observations we can deduce a simplified 
evolution kernel which reproduces with a good accuracy the exact one. Namely, the simplified 
anomalous dimension matrix reads 

ld J ) = ^ | 3 ^' + 3 ) " WO) " (— ^} (19) 



Mi: 


-> J 






-> J 





fc-A cj+l (j-j + ihj (J-j + ih 

and its eigenvalues coincide with the ones of Eq. ( |I7|) with a high precision, see Fig. ^. This 
anomalous dimension matrix corresponds to the kernel 

K = -^-{k(xi, X2\x[, x' 2 ) + k(x 2 , —x 3 \x 2 , —x' 3 ) + k(xi, —x 3 \x\, —x' 3 ) 

+ A(xi,x 2 \x' 1 ,x 2 ) + A(x 2 , -x 3 \x 2 , -x' 3 )j - \(3q5(xi - x' 1 )8(x 3 - x 3 ), (20) 

with k defined in Eq. ([□]) and 

A( Xl , x 2 \x[, x' 2 ) = 2 XlX ' 2 +f lX2 e° lu (x 1 , -x 2 , x, - x\) + 2^pLQ* n (x u -x 2 ). (21) 

JL-[JL 2 X^dj 2 

As is well known the tree level conformal invariance of QCD [^8[] allows for the diagonalization 



of the pair-wise evolution kernels in subchannels by means of the Jacobi polynomials p^ b ' Ua ' [E9 



j 



[P3| and express the two-particle Hamiltonians as a function of corresponding Casimir operators. 
Hence, we find 

(22) 

where j is a conserved quantum number related to the conformal spin of the composite two- 
particle conformal operator Ojk = 4>b{id a + idb) k Pj Ub ' Ua ^ ((d a — db)/(d a + <%)) <f) a , composed from 
elementary fields of conformal weight v — d + s — 1 (= 2 for gluons). Namely, the eigenvalues of 
the pair-wise Casimir operator of the collinear conformal algebra su(l, 1), J ab , on these states are 



S 



to 



22 . 5 



20 - 



17.5 
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Figure 3: The exact energy spectrum (|Tj) versus the selected trajectories generated by the 
approximate formula (|i~9l). 

J a b(J a b — 1) where J a b = j + \{y a + v + 2). Therefore, we can substitute the original problem of 
the diagonalization of the kernel (PH|) 



K = ~y H ({x t }\{x>}) + ^5( Xl - x[)5(x 3 - 4), 
by a quantum mechanical problem of the diagonalization of the Hamiltonian, H — > TC, 

n = n + v, 

where 

T-(.0 = hi2 + h 23 + h 31 , V = Ul2 + ^23, 

with the pair-wise interaction 

h ab = 2^(j ab ) -2^(1), ' 



fab 



(23) 



(24) 



(25) 



(26) 



a6 



Here the operator J ab is defined as a formal solution of the equation J ab = J a b{Jab ~ 1) and 

-2. ~2 -2 -2 -2 

possesses the eigenvalues J a t,. Obviously, [H, J ]_ = where J = J 12 + J 23 + J 31 — 9/4 is the 
total Casimir operator of the gluonic system and, therefore, it is the integral of motion. The 
spectrum of Ho part corresponds to the anomalous dimensions of the operator T. 
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-2 

The explicit form of the operators J ab depends on a basis where their action is defined. For 
our purposes we find it convenient to use the space spanned by the elements 9 k = r ^+v+i) \—> HI- 



Then, it is easy to read off their form from the action of the generators of the collinear conformal 
group on the elementary field operators. Finally, we have for the step-up, the step-down and the 
grade operators the representation 

i* = ( „ + 1)J + J .», = J. = l(„ + 1) + J ». (27 , 

The quadratic Casimir operator is given by J = J 3 (J 3 — 1) — J + J~ and equals j 2 = ^(z/ 2 — 1) 
for a one-particle state. The two-particle Casimir which enters Eq. (EBp is 



•4 = -t K+n)/2 a a a 6 ^+ ( ^ + ^ )/2 + \{v h - u a )6 ab (d a + d b ) . (28) 
2.3 Isotropic Heisenberg magnet. 

-2 

In this section we show that the Hamiltonian 7i$ possesses, apart from J , an extra integral 
of motion which makes the former exactly integrable. We will find that it corresponds to the 
generalized homogeneous XXX chain with sites of equal spin s = — — — | . To demonstrate 



this we use the formalism of the Quantum Inverse Scattering Method |36| |37], [38], HH anc ^ consider 
the lattice with equal conformal weights on each site vi = v. To prove the integrability of the 
Schrodinger equation 

H % = £ %, (29) 

let us define a matrix R a ^ acting on the product of the vector space of the site V a and an auxiliary 
space Vb- It satisfies the Yang-Baxter equation |57], 



The solution to this equation for the case when the dimensions of the quantum and auxiliary 
spaces coincide and which is of the most interest for us is given by |41], |32], |33 

T(J ab - \)T(i> + 1 + A) 

It is defined up to an arbitrary c-number function f(i>, A). In the case when the auxiliary space is 
two-dimensional the i?-matrix gives the Lax operator 

L a (X) = R,(\-l)=Xl + a i Jl (32) 



2 

where J = (J 1 , J 2 , J 3 ) with J 1 = § (j~ - J+) and J 2 = f (j" + J+). 



10 



Following a standard procedure we define the auxiliary and the fundamental monodromy ma- 
trices 

Ti(A) = L ai (A)L Q2 (A)L a3 (A), T b (X) = R ai , b (\)R a2:b (\)R a3 , b (\), (33) 

satisfying the .RTT-equation which is Eq. ([3(]) but with two i?-matrices replaced by the mon- 
odromies. It means that R(X—fi) intertwines the two possible co-multiplications of the T-matrices. 
This immediately leads to the commutation relations 



between the transfer matrices, 



MAi),t& 2 (A 2 )]_ = 



4(A) = tr b T b (A), 



(34) 



(35) 



which act on the total space of the spin chain <8>f =1 V ar 



The fundamental transfer matrix gives us the total Hamiltonian of the chain [^0], [IT], [42 

lntfr(A) h aia , 2 -\- h a2CL3 -\- h a3 ,a\i 



d 



dX 



(36) 



A=0 



with the two-site Hamiltonians 



h a ,b = R a M R a, b (°) = 2 ^ (J*b) ~ 2^(1), 



(37) 



coinciding with Eq. ( f2"6|) for the one-particle conformal weight v = 2 and provided f(u, A) 
r(i+A) r(^+i-A) • ^ ne ex P ans i° n °f the auxiliary transfer matrix in the rapidity A gives 



ii (A) = 2A + A J -3J 



where the total conformal spin and the 'hidden' charge are 



(3* 



r [j 1 + j 2 + j 3 y, Q = 2t l3k j[JiJl 



(39) 

respectively. Note that Q can be represented as a commutator of the two-particle Casimir opera- 
tors ( p8|) J ab = (J a + J b ) making use of the relation e^fc J\ J 2 jf = |[«/ 12) J 23]- which allows to 
rewrite it equivalently as [J 12 , J 23 ]_ = [J 13 , J 12 ]_ = [J 23 , J 13 ]_. 

For the later convenience we introduce a convention for the eigenvalues of the quadratic Casimir 
operator of the chain, namely, 

/7 2 =(j+^+l))(/+^(^ + l)-l). (40) 
For gluons at large J we have 77 = J + 4 + (9(J _1 ). 
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2.4 Conformal basis and recursion relation. 



The eigenfunction of the three-particle system can be found making use of the existence of the 
additional non-trivial integral of motion 



o- 



(41) 



where q are the eigenvalues of Q. 

The operator product expansion suggests that we must be primarily interested in polynomial 
solutions of the above equation. We can make further profit of conformal invariance of the system 
by looking the solution to this equation in the form of expansion w.r.t. a three-particle basis: 



where the basis vectors 

Pj-Ao 1 ,e 3 \6 2 ) = Q- j 



- (2, + 5) T(J + j + 6) px 



(j + 3) i+3 rV2( 2 J + 8) 



2j + 6 



(42) 



with 9 = ^-, (43) 

#23 



diagonalize simultaneously ^T|, ^] the total J and «7 13 Casimir operators 



a2 



JuVj-A !, W = U + 2) (j + 3) Vj^Oi, o 3 \e 2 ), 



(44) 
(45) 



and are normalized to unity (Pj'-j'(9 1 ,9 2 \0 3 )\Vj ] j(9i,9 2 \0 3 )) = Sj'j5j'j w.r.t. SU(1, 1) invariant 
scalar product 



<X2(01,02,0 3 )|Xl(01,02,03)> 



' d9fd9f , , 
' [l-OfBpY 1 - 1 
2m 



n 



(l-etOeY'-' x2(9i,e 2 ,e 3 ) X i(e 1 ,9 2 ,9 3 ). (46) 



\e e \<i 



The Hamiltonian H,q of the system exhibits cyclic permutation symmetry of the sites of the 
chain. Therefore, its eigenf unctions could change at most their phase under this transformation 



P%(9i,9 3 \e 2 ) = !Pb(02,0i|0s) = e l ^ (0i,9 3 \9 2 ) 



Then 



9^^ l f> 1 



(47) 



(4* 



for the "anharmonic" ratio of the conformal basis. Since the triple action of P leads to the system 
at its initial state, the eigenvalues of the operator are cubic roots of unity 



2 4 

<P = 0, r , -7T. 



(49) 
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Since this phase appears as a result of the cyclic shift of the chain sites it has an obvious inter- 
pretation of the momentum of the lattice. Obviously, it is the function of the conserved charges 
and is expressed through the equality 

E/ =0 HH2j + 5)7} 
Ey= ^(2j + 5)1} 

with Tj related to the expansion coefficients \Pj via Eq. fl5"2"|). It follows from this definition that 
= -¥>(?)■ 

The cyclic permutation transforms different three-particle bases and since each of them cor- 
responds just to a different quantum mechanical addition of three spins they are related by the 
conventional 6j-symbols as 

Vj-MMh) -j2 w AJ)'Pjd^W. (si) 

k=0 

The Racah coefficients fulfill the following properties: J2e=o WjtWkt = 5jk which is a consequence of 
the completeness condition £/=o \Vj-j){Vj-j\ = 1; T,e,m=o W j£ We m W mk = 5 jk coming from P 3 = 1; 
and VFj-fc = (-l) i+ *W fci reflecting the relation P 2 = P 13 PP 13 . 
Introducing new conventions 

= HYq^, (52) 



with 

s [ (j + l)4(J + J + 6) 2 (J-j + l) 2 ^' 2 
we can deduce from Eq. fl4T|) the recursion relation for the expansion coefficients Tj 



(53) 



gl} = ajT j+1 + bjTj-Xi (54) 

where 

(j + ]0(j_+_3) / 7 | ' i fc \ { 7 ' I 0\ I, (j +2)(j +4) ' i / 7 

^ = ^(2H^r (J+j+6)(J " j+2) ' 6j = ^(2 J T^r (J+j+7)(J " j+1) ' (55) 

and satisfying the boundary conditions 

T-i = T J+1 = 0. (56) 

Obviously, the coefficients Tj are related for positive and negative values of the eigenvalues of 
integral of motion due to P^QP^ 1 = — Q as 

U-q) = (-1)^). (57) 

And since PxsVj-jiOt, 9 3 \9 2 ) = (-l) J 'P J;i (^i, 6 3 \6 2 ) we conclude from here that P 13 & (q) = %(-q). 
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10 20 30 40 50 

J 



Figure 4: The log of the absolute value of the numerically diagonalized conserved charge q with 
q = trajectory being removed. The two possibilities of the description are clearly seen: from 
the top by the trajectories which behave at large J as q oc J 3 ; and from the bottom with the 
behaviour q oc J 2 . 

2.5 Conserved charge. 

The exact analytical solution based on the Eq. ([54]) is difficult. We, therefore, perform our analysis 
for large conformal spin J. 

There exists, however, the exceptional point q = when the exact solution can be easily found 
since Eq. Qol| ) is reduced to a two-term relation. Due to Eq. Qo"7| ) we have T 2 k+i(0) = with 
k — 0, 1, . . . , J/2 and from the condition of fulfilling the boundary condition 7j + i = we conclude 
that this solution exists only for even J. For even j = 2k we get 

r 2fc (o) - (-i) r (o) n — - ^ (^-^(^) fc To(0) - (58) 

Now we proceed with the rest of the spectrum. From the condition q 2 < ^ajbj of the existence 



of solutions to the recursion relation (|54"D we can easily deduce the critical values of the conserved 



charge for asymptotical conformal spins J — ► oo to be 



T3 t3 

< q < -V- (59) 



3\/3 ~ ~ 3^' 



14 



This equation gives the upper boundary for the quantized values of integral of motion shown in 
Fig. |j. This picture demonstrates as well that there exist two possibilities of parametrization of 
the spectrum by different sets of trajectories. 

The trajectories following the q = behave at J —>■ oo as q oc J 2 . Introducing the convention 
q* = q/rf and the continuum variable r = j j J the recursion relation can be reduced in leading- J 
approximation to a simple first order differential equation the solution to which reads 

1/2 /l_ r 2\< 



T 3 = Ji^(r), 



(r) = e-**' |r (3/2 + itf 



r 2 (l-r 2 ) 



(60) 



It is valid for whole interval of j namely for j 3> 1, J — j ^> 1 except of the vicinities of the 
reflection points j en d ~ 1, J where oscillating WKB solution is not applicable and we have to solve 
Eq. (53) close to the end points exactly. We have for J ^> J — j ~ 1 



2q*Tj = (J-J + 2)T j+1 + (J-j+ 1)2}-!, 



while for J ^> j ~ 1 



U + 1)0' + 3) 



T 



(j + 2)(j+4) 



2j + 5 ' 2j + 5 

To solve the first limiting recursion relation we make the following ansatz 



(61) 



(62) 



(63) 



Substituting it into Eq. (|6lD we get the equation 

d 



with the solution 



2A(1 - X)—v e (X) = (£ + l)u/+i(A) -(£ + 2)u/_i(A), 
dX 



vz{X) = C[X{l-X)]- 1 ' 2 {l-2Xy +2 . 



(64) 



(65) 



The solution to the second recurrent relation can be found by exploiting the cyclic symmetry 
of the problem and the known WKB solution (0). Using the Racah coefficients, Wjk, which relate 
the different three-particle bases fl5T|), we can write 



$ j = e i *'EW jk (J)$ k . 



(66) 



k=0 



From the recursion relation which is obeyed by the coefficients W (see e.g. Refs. [fJU], p2|) we can 



obtain the formula for the limit r 
1 



fixed, j = fixed, J — > oo, namely, 



W. 



jk 



77 



-l) J Wj(r), with Wj(r) = \p2 



r 5 (l - t 
AT" 



2\2 



1/2 



Cf /2 (2r 2 - 1) 



(67) 
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with the normalization coefficient Nj 1 = (3 ■ 4) 2 (2j + 5)/(j + 1)4. Then the solution for j ~ 1 
region is given by 

1 /o 1 

0(r)tWj(r) 



7) = z J '(-l) J e iv V2 



[0' + l)4 _ 


1/2 1 


_ 2 j + 5 _ 


/ 





T 3( 1 _ T 2)2]l/2- 



'Pi ? 



3 3 
2' 2 



(6* 



(69) 



(70) 



(71) 



Matching the WKB with the end-point solutions in the region of their overlap we get the 
quantization conditions for the charge 



A simple calculation gives us finally the result 

r. = ( i) J -j SL±lk 

3 \ ) ^ 

in terms of Askey- Wilson polynomials^ 

-j, j + 2a + 2(3 — 1, a — ix 



a + (3, 2a 



1 , 



Pj (x\a,/3) = i J 3 F 2 
which are orthogonal on the interval — 00 < x < 00 

dxpj (x\a,P)pk (x\a,P) \T(a + ix)T({3 + ix)\ 2 = 0, j ^ k. 



q* In i] = arg r ( - + iq*^ + ^ (2m + aj) 

where <jj = [1 — (— l) J ]/2. The first iteration at large 77 is 

. 7r 2m + oj 



(72) 



(73) 



6 In??- -0(3/2) 

with accuracy 0(\n~ 3 rj). 

Let us turn to the second way of describing the spectrum of q. The maximum of the eigenfunc- 



tions describing the trajectories which behave as q — ► J /v 27 is achieved for j max = ^ and in the 
vicinity of this point the system behaves as a classical one. To find the quasiclassical corrections 
we expand the recursion relation around this point j = ^= + Xy/j), and look for the solution 
to Eq. ( |5lD in the form of the series w.r.t. the inverse powers of J 

J 3 



*(A) =£<^)( A ) J ~ £/2 > and g(J,n) = -^E^W-T 
where 7} = <?(A). This leads to the infinite sequence of differential equations 

Z>(i)£( )(A) = 0, 
P (1) (1) (A)+P (2) (O) (A) = O, 

^(1)^(2) (A) + P( 2 )^(l)(A) + P( 3 )<P(0)(A) = 0, 



(74) 



(75) 



See also Refs. 



in connection to solution of the two-reggeon Baxter equation. 



16 



supplied with the boundary conditions <5(±oo) = and the differential operators given by the 
formulae 



The solutions ^(£)(A) = <j>(t){X) exp(— A 2 /2) to these equations are expressed by a linear combina- 
tion of the Hermite polynomials 

(O ) (A) = H n (A) , (76) 
0(D (A) = ^-n(n-l)(n-2)H n ^(X) + j(n-U+10V3)H n ^(X) 

io 4 
- l(n-17+10V3)H n+1 (X)-^H n+3 (X), (77) 

where n = 0, 1, . . . and gives the number of nodes of the solutions in the classically allowed region. 
Deriving these equations we have implicitly assumed n <C J in order the expansion (|74j ) for the 
charge to be meaningful. 

The quantization for the charge q immediately follows from boundary conditions. In this way 
we derive a few next-to-leading corrections to the charge 

s (0) H = i 

q {1) {n) = 7;~ n > 

q {2 \n) = -g--y™ + 3 n > ( 78 ) 
• • • > 

which compare quite well with the eigenvalues deduced from the numerical diagonalization of Eq. 
(H, see Fig. |. 

2.6 Eigenvalues of Hamiltonian Ho. 

Having found the expansion coefficients Tj, which completely specify the eigenfunction we can 
express the three-particle energy £ in their terms by exploiting the permutation symmetry of the 
system 

Wg )=2Re £ ^7^ + 5) f 7 > + 3, (79) 
1 q) E/=oHH2j + 5)7} 1 ; 
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Figure 5: The numerical quantized values of the charge (crosses) and WKB approximation 

mm- 



with the two-particle energy 



e(j) = 2^(j + 3)-2^(l). 



JO) 



Conventionally, taking the average of the Hamiltonian 7i w.r.t. the eigenfunctions \Pq we find an 
equivalent expression for the energy 



S (J,q) = 3 



S/ =0 e(j)|^ 
E/ =0 I^I 2 



These results provide the exact energy of the three-gluon system with the same helicity of all 
particles. 

Substituting the eigenfunction (|58"D corresponding to the special case when q = we easily 
obtain the lowest energy trajectory for even conformal spin J 

J \ „ , (J 



SoU 0) = 2^ (- + 3) + 2^ ( 2 + 2 ) - 4 ^(!) + 4 - 



82 



Next, using the WKB solutions found in the preceding section we can get the explicit result for 
the set of trajectories describing the spectrum from below using Eq. fl60D as 

'3 



S (J, q) = 4 In 77 - 6^(1) + 2Re tp (- - iq*^j 



(83) 
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J 

Figure 6: The numerical eigenvalues of the Hamiltonian 7i versus the analytical energy given 
in Eq. fl34|) for m — 0, 1 (odd J) and Eq. (|83"D for higher trajectories m = 3, 5 (even J), m — 7,13 
(odd J). 

where q* = g/77 2 and is valid with (9(?7 _1 ) accuracy. For the few lowest trajectories following the 
q = one we can write from here the approximate formula 

£o « 5 ) = 41n 2 - 4,(1) + 4 - (3/2) £ ( , n ,^; ( g 3; 2) ) 2 . (84) 

where we have used the explicit form of the quantized q from Eq. (^). For comparison with the 
numerical diagonalization see Fig. |6|. 

For the trajectories starting from the top of the energy spectrum we get using Eqs. (|7^j77| ) 
and (gl]) 

S(J,q) = 21ng-6^(l) + £>(J- 2 ), (85) 

where we have to substitute the explicit WKB expansion ( |78|) for the integral of motion. Evaluation 
of the higher WKB corrections to this equation is extremely complicated task because of the 
lengthy form of the WKB functions for i > 2 so that even 0(J~ 2 ) correction to energy 
becomes intractable. 



19 



2.7 Bet he ansatz and Baxter equation. 

As we have seen above the evaluation of the higher WKB correction to the energy fl85|) is a rather 
non-trivial task within the formalism used above as it requires cumbersome calculations due to a 
badly organized expansion in J: the actual expansion parameter for the eigenf unctions was \fl . In 
this section we will address this issue starting from the Baxter equation for the Q-operator. For the 
polynomial solutions the latter is equivalent to the usual algebraic Bethe ansatz [38|, [39|, [13| . 

To construct the diagonalized auxiliary transfer matrix ti (A) one defines a local vacuum state 
on each site \ue) so thatQ 

Ji\wi) = 0, Jt\wt) = \wt), (86) 

where from the definitions of the su(l, 1) generators we conclude that = const but it is not 
annihilated by the generator of the conformal spin J and thus is a legitimate choice. Then the 
vacuum of the chain is \Q) — <8>f =1 |u^) and the auxiliary monodromy matrix acts as 

T,(A)|^M« *«U )= ((^) 3 ° V (87) 
\ C(X) D(X) ) { * 0"^) ) 

where the * stands for a term whose explicit form does not matter for us. Hence, as we see the 
element B(X) is the annihilation operator while C(A) can treated as the creation one. Therefore, 
we can construct a special set of the states of the chain by acting with C(X) on the vacuum \f2), 
the so-called Bethe states 

\*j) = flC(\ t )\n), (88) 
1=1 

whose spin is Y^=iJfWj) — {J + \{y + [H . Making use of the exchange relations 

stemming from the RTT fundamental commutation relation one can easily deduce that these states 
are the eigenvectors of the auxiliary transfer matrix ti (A) = A(X) + -D(A) with the eigenvalues 



a+— n , \. + a -^h n r-r- » ( 89 ) 



m m m 



^ / £ =1 A — A^ V 2 / £=1 A — A £ 
provided the Bethe roots {A^|£ = 1, . . . , J} satisfy the Bethe Ansatz equation 

X m + n ( x m - A, + 1) = (x m - n (Am - A, - 1) . (90) 



8 Note that we can equally treat J7 as the creation and as the annihilation operators, respectively, as is 
done conventionally Q. This corresponds to the second choice of the pseudovacuum. In this case the local 
vacuum is defined as \ue) = Oj"^ 1 and the eigenvalue of Jf just changes the sign as compared to the Eqs. (|86|): 
Jf\u)t} — —^-\uje). In this case the Bethe states are = Y\.i—iB{\t)\Q) an d the eigenvalues of the auxiliary 
transfer matrix as well as the Bethe Ansatz equation remains intact. 
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On this condition the Bethe states become degree- J homogeneous translation invariant polyno- 
mials in Qg_. 



By virtue of Eq. (34) one can assume that the Bethe states are also the eigenstates of the 
fundamental transfer matrix. It is really the case and for the Hamiltonian of the lattice we have 
the expression in terms of Bethe roots W\ 53 1 



d 



(91) 



A=0 



Introducing finally the function 



Q(A) = J] (A - Xe) 



we have from Eq. (|89|) the Baxter equation |45], [46], [47 

v + V- 3 



ti(A)Q(A)= (A + 



Q(A + 1)+ A 



v + l 



Q(A-l) 



(92) 



(93) 



2 J " v 'V 2 
for the eigenvalues of the Q-operator where the ti(X) stands for the eigenvalues of the transfer 
matrix, rather then the operator, and reads 

ti (A) = 2A 3 + (V - 3j 2 ) A + iq. (94) 

The usefulness of the Baxter operator is that the energy of the system can be rewritten concisely 



in its terms making use of the known Bethe ansatz representation |24 

So 



Q> to) Q> (-^ 



+ 6^(z/ + l) -6^(1), 



(95) 



provided f(v,X) is given by the expression after Eq. (|37D . Thus, the knowledge of the Baxter 
operator as a function of the conserved charges allows to find immediately the energy as a function 
of the latter. The Eq. (P3J) will be the main object of our analysis in this section. 



2.8 Quasiclassical expansion. 

Let us analyze the Baxter equation for large values of the total conformal spin, J, of the three- 
particle system as was done before when we have dealt with the recursion relation (|55J). In these 
case the system behaves as a quasiclassical one, and we can apply appropriate methods used 
previously for the analysis of the Toda |48| and s — — 1 Heisenberg spin chain [^6 . 



To get rid of the large factor J in the transfer matrix let us rescale the spectral parameter as 
A —>■ JX and introduce instead of Q and ti the functions 



£(A) = Q(JA), T(X) = (JX)-H h (JX) 



(96) 
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Then Eq. (p3|) takes the form 



A + ^±lj-^\(X + J- 1 ) + (a - ^±1 j-^ 3 fa - J' 1 ) = A 3 T(A)£(A), (97) 
where T{iX) is a real function. Moreover we introduce conventionally the WKB function as 

£(A) = e JS( ^ x \ (98) 

so that the energy is expressed now as 

£o = S' {^J- 1 ) - S' + ^{v + 1) - 6^(1), (99) 

and, thus, for the limit J — > oo we can primarily be interested in the small A asymptotics of Eq. 
fl9~T|). We look the solution to this equation in terms of series w.r.t. inverse powers of the conformal 
spin 

oo oo 

S{\) = Y / S w (\)J- e with T(A) = ^%(A)J- £ , (100) 

1=0 1=0 

where 

^ 0) (A) = 2 + A- 2 + ^g (0) A- 3 , ^)(A) = 8A- 2 + -^=g (1) A -3 , 
^)(A) = yA- 2 + -^=g( 2 )A- 3 , ^(A) = -^>A" 3 , for t > 3. (101) 

In leading order approximation we have the equation for SL^ 

2coshS| 0) (A) = 7 w (A), (102) 

with the solution 

S{ 0) W = 21n|^7(o)(A) + i + ^^(A) - || , (103) 
which has a correct asymptotics for A — > oo since 5(A) A ~°° In A as follows from Eq. fl92|) because 



Q(A) ~ A . In the small A region we get from Eq. ( Jl 03 ) 



^(o)(A-±0) = ±ln-^g(°)A- 3 . (104) 

From here we can immediately find that for large conformal spin Sq{J — > oo) = 6 In J and coincides 
with the result flS5f ). Unfortunately, due to the fact that for the calculation of the energy from 
Eq. ( p9|) the function 5* enters in the point ~ J -1 in order to estimate a non- leading correction 
we have to solve an infinite series of differential equations stemming from Eq. (0). This is, of 
course, not feasible. Therefore, we have to perform an effective resummation of these series |26|. 
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To do this we notice that from Eq. ( |104| ) is follows that for ReA > we have exp ^5| ^(A)^ ^> 
exp (— S'^(X)^ and one can easily convince oneself |2(J that up to 0(X 6 ) the small A expansion of 
Eq. (0), i.e. Eqs. (|114j) , can be substituted by the simplified equation 

^ X + J ^-X!(x + »-±lj-iY 3 T { X). (105) 



«A) V 2 

Similarly for ReA < we can omit the first term on the l.h.s. of Eq. (]97D . Now these simplified 
equation can be solved exactly as 

S'(X) = 3^( JA) + 3^ ( JX + ^r^) + 7lnT(A), for ReA > 0. (106) 

V 2 / exp(J _1 a) — 1 

And analogously for the ReA < obtained from this equation by J — > — J and changing the sign 
of the third term. One can readily recognize the differential operator as the generating function 
t/(e t — 1) = J2eLoBgt e /i\ for the Bernoulli numbers of the first order [j5(|. 
Finally, we have for the energy 



£ = QiP (^-) -6V(l) + 2Re 



lnT(A) + 0(J~ b ). (107) 



exp(J -1 <9) - 

Making use of the representation of the auxiliary transfer matrix in terms of its zeros 

3 

ti(\) = 2l[{\-r e ), (108) 



one can easily obtain a concise expression of the energy (cf. |26| 

£ = 21n2 + 2Re]T(V(^-r,) - V(l) ) + 0(J~ 6 ). (109) 



f — i ^ 

2.9 Quantization of q and energy. 

Expressing the energy directly in terms of the integrals of motion we have to find their quantized 



values. Although it was done in section |2.4| from the study of the recursion relation we address 
this question here from the point of view of the Baxter equation. From Eq. (p2|) we get for S'(X) 



S'(X) = ^(JA - A*)" 1 , (110) 
e=i 

and thus 

^(A) = l, (111) 

c 

for a contour in the complex A-plane which encircles all the roots Xi/J. From the leading WKB 
solution we conclude that the zeros of Q(X) are accumulated on the interval |T(A)| > 2 on the 
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imaginary axis. For the case at hand we have two intervals of instability for the classical motion 
|49| f4"8[ , |2~6| and they are determined by the system of cuts [Ai, A 2 ] U [A3, A 4 ]. Therefore, we have 
n zeros on [Ai, A2] and J — n on [A3, A4] and, thus, 

f^l sw = r (112) 

Ci 

with the contour C\ which encircles the interval [Ai, A2] (see Fig. |7|). It is obvious that < n < 
J. Let us assume for the following that n ~ 0(J°) and substitute the expansion (|100|) into 
quantization condition (|112|) to get the sequence 

f£iW=°. f^S' m W = n, = 0, for I > 2. (113) 

Ci Ci Ci 

The functions 5*^ can be easily deduced from the Baxter equation (|57|) and read 
/ ^(i) 1 // / 9 ! 

1 /// 1 / // 9 i // 



and S'^ has been found before in Eq. (|103|) . From the quantization condition for SL-s and Fig. 
it is seen that the function S'^ should not have singularities on the interval [Aj,A2] which is 
the case provided the latter shrinks into the point A cr j t . These critical values are defined by the 
equations T/ Q JX cr i t ) = and 7(o)(A cr it) = —2 with the solution 

1 i 
g C rit = q {0) = -, A crit = —^7g ( 115 ) 

where we have restricted ourselves to > 0. 

To proceed further we consider an infinitesimal interval around A cr it 

A 2 -Ai=ie, Ai = A crit - ie/2, A 2 = A crit + ie/2, (116) 

to ensure the presence of singularities of S'm\ on it and to fulfill the second quantization condition 
( |113|) . Here 7( )(Ai) = T(q^(X 2 ) = —2 as is seen from the Fig. |7|. Expansion of the rescaled transfer 



matrix around the critical point A = A cr i t + \e{2x — 1) gives 

T {0) = -2 + U4x(x -l)e 2 + 0(e 3 ). (117) 
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Figure 7: Typical form of the rescaled transfer matrix 
Inserting it into Eqs. ( |103| , |114| ) we get 



for < < i. 



ITT ■ 

i 



- 12iy/x(x- \)e + 0{e 2 ), 
7(i)(A cr i t ) i 2x-l 



2As 



X[X 



As x(x — 1; 



'118) 



Substituting the last equation into the quantization condition is § Ci = n we can deform 

the integration contour away from the cut [Ai, A2] reverting its direction and take the residue at 
infinity x = 00 to get 

- t£ Res — \—T { 1) (A crit ) + -| =n. (119) 
x =°° xe 1 24 w 2 J 

Using explicit form of from Eq. ( |101| ) one can immediately convince oneself that this is exactly 
in Eq. (|78]). 

Since the rescaled transfer matrix possesses the singularities only on the interval x G [0, 1] 
and infinity, we can always deform the integration contour and evaluate the residue at x = 00. 
Thus, for evaluation of the quantized q one has to keep only (xe) _1 -terms in the expansions. 
This simplifies considerably the task of calculation of the higher WKB corrections (and is easy to 
formalize) as compared to the method based on the solution of the recursion relation in section 
For instance, to find q^ we have just to keep more terms in the expansion in Eqs. ( p.l7|J118| ): 



(0) 



4 + 36(x £ ) 2 + 192v / 3(xe) 3 + 2160(xe) 4 + 
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^T a) = (g« -4)+2V3 (3qU - 8 ) (are) + 72 (g« - 2) (xe) 2 + ..., (120) 

which was done having in mind that we will be interested only in the residue at x = oo, xe = fixed. 
And, thus, we have 



S' {Q) = i {vr + 12(xe) + 32V3(xe) 2 + 304(xe) 3 + ...}, 

2,0-7 , 10,«-39 + 2 (669(1) _ 239)(i£)+ 1 (m) 



{ 2{xe) y/3 
to the same accuracy. Therefore, we have finally the necessary terms 



S'm = ( ^)(A crit ) - \ - %<P + Si + • • • • (122) 



^-~e\2i I{2)KAcnt/ 3 V ; 3^ ' 36 J 

From the requirement that this residue must equal zero we get the result in Eq. (|78|). 

One can generalize in a straightforward way this consideration for q^ with £ > 3, however, we 
will not pursue this goal here since in the consequent consideration of the perturbed Hamiltonian 
(Hp) we will be able to calculate the energy only up to 0(J~ 3 ) accuracy. The main motivation 
to study the Baxter equation was to express the energy directly via the integral of motion which 



allows to find the former in a very economic way. Applying Eq. ( |107| ) we easily get 

(o)j3 OoW 2 I 4 (V 1)N ) 2 - ll\ 

£ » - 2,n -vr + + T>[h - j y^J - 6 *' (1) + ° (J »• (123) 

where we have to plug in the explicit expressions for the expansion coefficients of the charge (|78|). 
2.10 Perturbed Hamiltonian. 

Let us turn to the perturbed Hamiltonian (p4|) whose eigenvalues are responsible for the scale de- 
pendence of the three-gluon correlator which contributes to the transverse spin structure function 
g2- Of course, the problem now is more complicated since the interaction is not integrable. But it 
is easy to figure out that the non-integrable addendum V is small for a bulk of the spectrum and 
this allows to consider it SIS db perturbation. 

To start with, let us introduce the more general Hamiltonian 

H a = H + aV (124) 

specified by a coupling constant a. Obviously H a= o = TCq and 7^0=1 = H. Considering the flow of 
the energy levels given in Fig. ^| we conclude that the upper part of the spectrum is not affected 
significantly by the perturbation V and it represents a small correction. At the same time it 
modifies in a sizable way the level spacing for the lowest trajectories. 
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Figure 8: Flow of the eigenvalues of the 7i a for J = 35. The horizontal dashed line shows the 
position of the level £o(0). 

The present problem does not manifest the cyclic symmetry but only w.r.t. the permutation 
of the first and third gluons. Therefore, due to antisymmetry of the three-gluon operator we will 
be interested in the eigenfunctions, \&(~\ of the Hamiltonian (|24|) 

(125) 



which are odd functions w.r.t. the permutation of the end gluons 



(126) 



For the top of the spectrum we have |V| <C \7io\ for large J and, therefore, as a first approximation 
we take 

with #l j - ) (q) = -^{%(q)-%{-q)}. (127) 
Then the matrix element of V in this basis reads 



V, 



g ,q 



(^ ( " ) (g / )l^ ) (? , ))(^ ) (9)l^o H (9))l -1/2 (^"(gOlvi^Cg)) 



'1281 



with the explicit eigenfunctions being plugged in there, it gives 



-1/2 



V^ = -8(£|^(g0rEl^)P 

Vi=o i=o 
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x E 7 ■ I ^ L {cos (^') - <?(?)) " cos ( V (g0 + <p(q))} . (129) 

U + 2)0 + 3) i- 

It defines the correction to the eigenfunctions 

where q = q(n) and q' = q(n'), and correction to the energy which can be easily evaluated making 
use of Eqs. ( f76[ , f7?|) with the result 

24 48 

V q , q = - T2 - T3 (2n-3) + 0(J-'). (130) 

Therefore, it affects the spacing of integrable Hamiltonian 7i only at order 0(J~ 2 ) while the net 
spacing is 5£ (q — > J 3 / '\/27) oc (9(J _1 ) and the former is a small effect for large J. 

As is seen from the picture of the flow of the energy eigenvalues of the Hamiltonian 7i(a) 
the dependence on a is linear. Thus, going to the limit a — > oo, we can consider the integrable 
Hamiltonian as a correction on the background of V. Calculating its matrix elements w.r.t. 
eigenfunctions of V and extrapolating the result to the point a = 1 we can estimate the spacing 
for the low trajectories of the 7i. 

The eigenfunctions of the V kernel which diagonalize it in the limit of large J are 

w v jI~ ] = -j={Vj im (e 1 ,e 2 \e 3 )-v J ? n (o 3 ,o 2 \e 1 )}, (131) 

with required antisymmetry w.r.t. the P 13 permutation. They are normalized as (^J.^ I^j.^ ) = 
8m,m' + 0(J~ 2 ). We have the eigenvalues 

V(m) = (^IVI^V^SV^) = - (m + 2) 4 (m + 3) + °( J ~ 2 )- ( 132 ) 

The 0(J- 2 ) terms come form SV(m) = V(m) + (m+2) 4 (m+3) = ^jg^ - 4Eio which 
can be estimated as follows. In the sum the main contributions come for the regions m, £ ~ 1 and 
< t < 1 with £ = rJ. For the latter we can replace the sum by the integral, J2 ~~ * J Joi an d 
use W from Eq. (|6T| ) so that one can immediately conclude that sum oc J -2 . For the small fixed 
m and £, and asymptotical J, as can be deduced from the recursion relation obeyed by W, the 
Racah coefficients are given by 

W jk = ^^(2A: + 5)(j + 2)(j + 3) + 0(J- 4 ). (133) 

and, therefore, we have finally for the matrix element SV(m) oc 0(J~ 2 ). This estimation is an 
example of the asymptotic addition law for anomalous dimensions of scalar multi-particles com- 
posite operators [^TJ, according to which, for the present case, the limiting point of the spectrum 
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at large J is given by the naive sum of the anomalous dimension of a field and of a two-particle 
composite operator. One can immediately recognize in Eq. (|132|) the anomalous dimensions of the 

%y 



two-particle local composite operators in ^jt-, -theory. 



Now the correction to the spectrum ( |132j ) is given by the a l (&j.^\Ho\ l l/j^) and reads 

(^noW^/i^Wj-^) = <rn) (1 - (-l) m W mm ) + E (2 " (-if) e(i)W m£ W me . (134) 
Making use of the explicit asymptotics of the Racah coefficients we get finally (up to 0{ri~ 1 )) 
£ a (m) = aV(m) + 4 lnr? + 4^(m + 3) + 2ip(m + 5) - 2^(2m + 5) - 2^(2m + 6) - 6^(1). (135) 

Since the dependence on a observed for the energy flow is linear we can extrapolate this result 
to the point a = 1 and estimate the spacing for the lowest trajectories of £. For instance for 
m = we get £i(0) — £o(J, 0) = —0.46 which is rather close to the numerical value found from 
the explicit diagonalization. 



Making use of these results we design the following formulae which describe well (see Figs. [10 
^) the spectrum of the perturbed Hamiltonian 

£(J,m) = £ (J,0) - A(m), for m = 0, 1, (136) 
£{J,m) = £ (J,q(m)) - A(ra), for m > 2. (137) 

Here £o(J, 0) is the ground state energy (|^) and £o(J, q(jn)) is Eq. (|^) with the q{m) trajectories 
deduced from the quantization condition (^). The function A is 

A(0) = 0.54, A(l)=0.08, A(m > 2) = (S - 5(m - 2)) 9 (S - 6(m - 2)) , (138) 

with #o = 0.15 and § = 0.01. It serves to generate a shift of the non-perturbed energy trajectories. 
The top of the spectrum coincides with the £q trajectories. 



3 Quark-gluon sector. 

Although we have started with the pure gluonic sector, chronologically, it is the non-singlet quark- 
gluon-quark sector which has been first addressed within the context of the transversely polarized 
structure function ^(^Bj)- The first study which correctly identified a complete basis of operators 
was by Shuryak and Vainshtein and they have calculated a lowest anomalous dimension [ j52]| . 
The evolution equation for the correlation function Y in Eq. (0) which is defined as the Fourier 
transform 

Y( Xl ,x 3 ) = f ^l^e^-^ihiy^O,^) (139) 
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Figure 9: The exact energy eigenvalues versus the approximate formulae ( |136| ). 
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Figure 10: Same as in Fig. |9| but for large J. The spectrum is cut at the top where it coincides 
with the eigenvalues of the integrable Hamiltonian 7i - 
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of the C-even combination 

y(ni, 0, k 3 ) = - 0, k 3 ) + ~S(-k 3 , 0, -«i)} , 

of the three-particle quark-gluon-quark operators 

^(ki, k 2 , k 3 ) = -^(K 3 n)ig^ + \iG^ + (K 2 n) ± 7 5 G^ + (/t 2 n)] V»(«in), 



(140) 



(141) 



has been derived in Refs. || |53| as well as the anomalous dimensions of the local operators were 
calculated. Consequently this was redone in a number of studies [^7], |54| with almost the same 
results. 



3.1 Evolution equation. 



In spite of the available results |k| let us discuss the evolution equation from our present point 
of view. The function Y can be decomposed according to this discussion as Y = + S + ~S. The 
operators <S are related to each other by the charge conjugation and, therefore, they evolve 
autonomously. The evolution equation for the functions S has the same form as Eq. (|6]) . For + S 
the total kernel is a sum of the ones in subchannels 



K 



+s 



qq Kf 8) ( Xl , -x 3 \x[, -4) + m K[ 3) ( Xl , x 2 \x[, x' 2 ) + q9 Kl 3) (x 2 , -x 3 \x' 2 , -x' 3 
^PoS(xi - x' l )8{x 3 - x' 3 ), 



(142) 



and their explicit form can be found in Ref. [21]. Analogous expression holds for ~S with T and 
V subscripts being interchanged in the quark-gluon kernels. Since these kernels are equivalent we 
consider in the following only + S. The pair-wise kernels contain both Cp and Cp — Ca/2 colour 
structures and this presents one of the sources of complication to find the analytical solution. 
Since for QCD (Ca — 2Cf)/N c = 1/9 is small we can anticipate that large N c approximation 
presents a good (within 10% accuracy) starting point for solving the three-particle problem. 

Taking the Mellin moments (|T3J) we can write the evolution equation in large iV c limit for the 
moments 



d 



dlnQ 2 



Ait 



(143) 



fc=0 



with the large- N c anomalous dimension matrix [53 



1% = 8 jk (^(j + l) + ^(j+4)+^(J-J + 2) + ^(J-j + 3)-4^(l)-^ 

(k + 2) 2 (J_A; + 1) 2 

— f, 



k,j+i 



(k-j)(J-j + l) 2 



(144) 
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For large N c we have "K%\ ocl/N c and the eigenvalues of the other kernels read, up to terms 
suppressed in N c , 

f dx 1 dx 2 Pj 1 ' 2) ( Xl X2 ) { J® 1 (x 1 ,x 2 \x[,x' 2 ) 
J \xi + x 2 J I wKL, I 



j\/e r vo' + 1) + v»(j + 4) - 2^(1) - § i p(1)2) 

2 \V0' + 2) + V(j' + 3)-2^(l)-|J J U + ^/ 



Therefore, similarly to the previous study we can replace the original problem of the diagonaliza- 

2 2 

WP = (146) 



tion of the three-particle kernel K +s = — > ^f-Ti by the solution of the Schrodinger equation 



with the Hamiltonian 



H = h 12 + h 23 - | (147) 



where the pair-wise Hamiltonians 

hl2 = V (3l2 + (jit ~ l) -2^(1), ^23 = V (^23 + 7j) + V> (^23 " ^) "2^(1), (148) 

can be easily read from Eq. (|145f) . 

3.2 Inhomogeneous open spin chain. 

It turns out that the problem (|146[ - |148|) is exactly solvable^]. Let us consider the inhomogeneous 
spin chain which is characterized by 'inhomogeneities' 5t. The monodromy matrix is defined as 
before but with the shifted spectral parameters 

21(A) = i? ailb (A - 5i)R a2 fi(X - 8 2 )R a3tb (X - 5 3 ). (149) 



The generating function of the conserved integrals of motion, i.e. the transfer matrix, reads |57] 

t b (X) = tr b K + (\)T b (\)K-(\)T b -\-\), (150) 

with boundary reflection matrices K which we set equal to unit matrix in order to ensure the 
conformal invariance of t b . The matrices (|1 50|) obey the commutation relation (0). The shifted 
and renormalized auxiliary transfer matrix nf=i {(A — <^)(A — 5e — 1) — j^jti (A — |) — > ti(X) 



U(X) = tri L ai (A - 5 1 )L a2 (X - 5 2 )L a3 (X - S^L^-X - 5 3 )L t 02 (-A - SJL^-X - 5 1 )a 2 , (151) 



3 See Ref. p5| for spin models appeared in the solution of the quark-gluon reggeon interaction. 
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which is an even function in A, has the following expansion in the rapidity^ 

fi(A) = 0(A) - (4A 2 - 1) (A 2 - 5\ - jl) f - i(4A 2 - 1)Q S , (152) 
where i7(A) is the c-number function 

12(A) = -2 f[ (A 2 + f t - 5f) + (4A 2 - 1) £ jg (A 2 + 8\ + «5f - *?) + (4A 2 - l)J 2 j 2 , 



i2 



J is the total Casimir operator and Qs is the non-trivial integral of motion 

Qs = [f 12 , J 23 ] + - 2(5\ - %)f w - 2(5 2 3 - S\)j\ 2 + mt l3k J[P 2 Jl (153) 



As will be clear later on, only 5 2 = spin chains appear in QCD within the twist-three context. 
Therefore, we set in what follows 5 2 = 0. 

The Hamiltonian which commutes with the charges found above is deduced from the funda- 



mental transfer matrix and reads |57j 7~Ls — ^12 + ^23, with the two-site Hamiltonians expressed 



in terms of the fundamental i?-matrices h aim = R aua2 (—S)R' ai a2 (—5) so that explicitly 

h 12 = V (Jn + St) + 1/) (J u - 5i) - 2^(1), h 23 = i, (j 23 + S 3 ) + i> (J 23 - S 3 ) - 2^(1), (154) 
provided we set in the bundle R a b(X) ( ]3T| ) the function f(u,X) as was done after Eq. fl3"7|). The 

~2 ~2 *-2 "2 

commutation rule for the h aiia2 with Qs is [Qs, hi 2 ]_ = 2[ J 23 , «7 12 ]- and ^23]- = 2[«7 12 , ^23]- • 
So that the property H]_ = is obvious. 

Setting <5i = 3/2, <5 3 = 1/2 we get the QCD result QH8|) and Q 5 charge of Ref. |9|. The 



conformal weights of the end quarks are V\ = u 3 — 1 and v 2 = 2 for the gluon. 

3.3 Master equation and integral of motion. 

As before we can solve now the simplified problem 

Qs* = qs$, (155) 

instead of Eq. ( |146| ), where Qs is Eq. ( |153| ) with <5's set earlier. The main steps of the solution of 
this problem are the same as we have pursued in the preceding section and in Ref. [^2] so that we 
just outline them briefly. 

~2 -2 

First it is convenient to take the polynomials, which diagonalize simultaneously J and «7 12 , 



Vj-jVi, m) - 7f (j + + 3) T(2j + 4) rV2( 2 J + 6) ^ 6 * Fl { 2j + 5 



9), (156) 



10 



We would like to thank G. Korchemsky for a discussion on this point (see also recent Ref. 
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with 9 = , as basis vectors since the off-diagonal elements of the Qs in this basis will be the 



same as for Qt in so that only diagonal elements will be slightly modified. Expanding the 
eigenf unctions \P in the way it was done in (|42|) 



(157) 



it is easy to show that Eq. ( |155| ) is equivalent to the three-term recursion relation 



(2j + 3)T i+1 + (2j + 5)7}.! + 6j (2j + 3)(2j + 5) ([Q s (|, J)]^ - 9s ) 7) = 



(158) 



with boundary conditions T_\ = Xj +1 = 0, and where 



3,3 



x 



j+ D( j ' + C l+u 



1) ( J + l) + l ( J + i)( J + 

2) l J+ 2j + 4( J + |)( J + |) 



(159) 



Here the new expansion coefficients are introduced as follows 

"(j + l) 3 (j + 3) 3 



#3 = ft T i' 



with 



63 



(J-j + l)(J + j + 5) 



-1/2 



(160) 



U + 2) 3 

However, due to the loss of any symmetry of the integral of motion w.r.t. the permutation of the 
sites it is instructive to study the recurrence relation based on the expansion of the eigenfunctions 
in the basis Vj-j(9 3 , 9 2 \9\) 

(161) 



3=0 



where 



T(2j + 5) rV2(2J + 6) 



2j 



Introducing 



with 6* = — — . 

9n 
(162) 



</'; - ^ / — / • with 9 = [(j + i)(j + 2)(j + 3)(J-j + l)(J + j + 5)r 1/2 , (163) 

we have the recursion relation for Ej 

(2j + 3) + (2j + 5)^_! + S 2 (2j + 3)(2j + 5) ([Q s (|, I)],-, - 9s) H, = 0. (164) 

From the condition of existence of a solution to the recursion relations which satisfies the 
boundary conditions we find that the allowed values of qg asymptotically lie in the band < 
qs/J 4 < §■ The upper boundary is achieved at j max = -j^J ■ In complete analogy with Qt [ f22|| , 
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the spectrum of Qs can be described by means of two different sets of trajectories. The first ones 
which behave as J 2 and another ones as J 4 , at large J. 



Let us turn to the first possibility. One can immediately find from Eq. ( 158 ) the exact lowest 
trajectory for the charge 

&"*W = (j+l)(j+l)+l- (165) 
For other levels we construct an effective WKB approximation. For the case at hand the classical 
motion is allowed for the whole interval of j except of the vicinities of the reflection points j e nd ~ 
1, J, which will be parametrized by the continuous parameter r = j/J. Introducing the new 
function as Tj = J(— l) J t;(r) and q~g = qs/J 2 , we have from (|158|) the differential equation 



r 2 (l - r>"(r) - r(l - t>'(t) + 2(g£ - l)v(r) = 



the solution to which reads 



v 



V) = C (+) r 1+2iris 2 F 1 2 



1 +ir)s, -§ + «7s 



1 + 2ir) S 



r 2 + 



/nr(-) 1+2m,s p ( \ ~ % VS, ~\ - IVS 

2 I 1-2% 



where rjs = \^J2q* s — 3. Analogously for Sj = J(— l) J ^(r) we have the equation 



r 2 (l - r 2 )£"(r) - r(l - T 2 )£'(r) + 2(g* - 1 + 4t 2 )£(t) = 0, 



with the solution 



f (r) = cWr 1+ai " s 3 Fi 2 



1 + 2ir/ 5 



r 2 +C H T W% 2fl 2 



1 - 2ir) S 



(166) 



(167) 



(168) 



(169) 



Here the arbitrary complex constants C^' and (7^ have to be fixed from a sewing procedure 
with the solutions \Pj and $j close to the boundaries. 
For J 3> J — j ~ 1 we can easily solve (|158[ ) 



170) 



To solve the recursion relation for J ^> j ~ 1 is more difficult but we accept the same strategy as 
for Eq. (|62j) . However, since the permutation symmetry is lost we have to modify the procedure 
accordingly. Namely, the expansion coefficients in Eqs. ( |157| , |161p are related to each other by 



j j 

*j = E w kj $ k , $3 = E w jk $k, 

fc=0 fc=0 



with the Racah coefficients W jk , P J;j (# 3 , 2 \0\) = Y. J k=QW jk (J)Vj- k {d l ,d 2 \d- i ), which can be read 
off from the general result of Ref. ||22|| . Then, for the region in question the solutions can be 
written as 



T 



(-l) J (i + 3) /drPf 2) (2r-l 



(j+2)v(t) 



(172) 
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Evaluating the integrals we obtain (see Ref. 



r^j-iyV + W + Q fe+> 

X4F3 1 " 



2*1 



1 + 2zr/ 5 
1 



1 



; - r(| + »^)r(| 



— I — J + ^5, I + J + «7s, 1 + 2i^5 



r(|-^ s )r(f +i + z?7s 

(173) 



and 



-• = (-l) J (j + 2) few 



2*1 



1 + 2ir/ 5 



1 +(-r 



,- r (i + ^) r ( 



x 4 F, 



lT]s 

I * ~ i 

r (I - r (f + j + zr/ 5 



+ 



<4-T3 I o t- 

Comparing Eqs. (|167 ) and ( |170|J173| , |174| ) in the region of their overlap we find 



-Vs] 



c(+) 


„E f 3-"7s,-!-«7s 


1) 


£(+) 




1) 


c(-) 


n P ( k+ ir >3,~j+iVS 




c(-) 


2- 1 1 V l+2i?7S 






(174) 



(175) 



and the quantization condition for the eigenvalues qs (cf. |jT9| 

C(-) C<+) T 2 (I + Z77s 



4r/5 In J — arg ■ 



2rcm, 



(176) 



C(+) CHpff-^) 
where m G Z + . 

For the upper part of the spectrum qs oc J 4 , the classical motion is concentrated in the vicinity 
of the point j max . Looking the solution to the recursion relation in the form of the series ([74]) we 
find immediately that the differential equations resulted from these are almost the same, to the 
required accuracy, as for the chiral-odd sector so that we can immediately adapt the results of 
Ref. f22|. Therefore, the WKB expansion of the charge reads 



Qs(J) 



00 m 



(177) 



with the expansion coefficients 

?s(n) = 
q { s\n) = 



1 

2" 
6 



V2 



y/2- 



n. 



(2)/ \ 
Qs ( n ) 



(403 - 72V2) + - (ll - 72V2) n + 



11 



-n 



:i78) 
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with n — 0, 1, It differs form the qr expression in Ref. JZIJ only at the level of 0(J 2 ) 

corrections. 



3.4 Energy spectrum. 

The knowledge of eigenf unctions allows to find the energy of the quark-gluon-quark system. For 
the trajectories starting from the top of the spectrum we have 

£(J, q s ) = \nq s /2 - #(1) - ~ + 0(J^ 1 ), (179) 

with qs taken from Eq. (|177| ). 

For the levels behaving as £ oc 2 In J we can find the exact lowest trajectory |58| corresponding 
to the charge ( |165| ) 

£{J) = ip(J + 3) + if>(J + 4) - 2^(1) - 1 (180) 
The remainder of the spectrum is described by the formula [|T^, ^] 

S(J } q s ) = 2 In J - #(1) + 2 Re V (J + zry s ) - | (181) 
and it compares quite well with the explicit numerical diagonalization of Eq. ( |144j) as shown in 



Figs. 11, 12 



4 Conclusion. 

In this paper we have demonstrated that the one-loop QCD evolution equations for the three- 
gluon and quark-gluon-quark correlation functions can be reduced to (almost) integrable one 
dimensional lattice models. The appearance of the integrable interaction deduced from the Yang- 



Baxter bundle (|3l|) , within the present context, is a consequence of the logarithmic behaviour of 
QCD in collinear regime since the Mellin transform w.r.t. momentum of the emitted gluon from 
a particle diverges as In J for large J, i.e. for the kinematics with negligible momentum transfer 
to a gluon. For finite J it stems from the ip( J)-function, i.e. the one which appears in local 
Hamiltonians. It is well known that in the region of momenta mentioned above, QCD can be 
replaces by the theory of eikonalized fields which should possesses, hence, more symmetries than 
the original QCD Lagrangian. 

Integrability allows to disentangle the energy spectra of the three-particle systems and an 
additional hidden charge serves to count the trajectories giving rise to different scale dependence 
of the corresponding multiplicatively renormalizable components of the correlators. The explicit 
analytical solutions have been found in the form of WKB type expansions w.r.t. the total conformal 
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Figure 11: The analytical trajectories from the Eqs. (180, 181| ) and the numerically diagonalized 



anomalous dimension matrix (144). 
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Figure 12: Same as in Fig. |TT] but for large J. The spectrum is cut from the top where it is 
almost continuous. 
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spin J of the system. A few first terms in the series gave a good accuracy as compared to the 
numerical results obtained by a brute force diagonalization of the anomalous dimension matrices. 
Moreover, our analytical expressions are valid up to a very low J and, therefore, this allows to apply 
the solutions for realistic calculations of the effects of scaling violations in physical observables. 

We would like to thank V.M. Braun and D. Miiller for a conversation at the preliminary stages 
of the work and especially G.P. Korchemsky for a number of discussions during the visit of the 
present author at LPTHE (Orsay) in May 1998. We acknowledge the support from the Alexander 
von Humboldt Foundation. 
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